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Abstract. A description of tlie event liorizon of a perturbed Scliwarzscliild 
black iiole is provided in terms of the intrinsic and extrinsic geometries of the 
null hypersurface. This description relies on a Gauss-Codazzi theory of null 
hypersurfaces embedded in spacetime, which extends the standard theory of 
spacelike and timelike hypersurfaces involving the first and second fundamental 
forms. We show that the intrinsic geometry of the event horizon is invariant 
under a reparameterization of the null generators, and that the extrinsic geometry 
depends on the parameterization. Stated differently, we show that while the 
extrinsic geometry depends on the choice of gauge, the intrinsic geometry is gauge 
invariant. We apply the formalism to solutions to the vacuum field equations 
that describe a tidally deformed black hole. In a first instance we consider a 
slowly-varying, quadrupolar tidal field imposed on the black hole, and in a second 
instance we examine the tide raised during a close parabolic encounter between 
the black hole and a small orbiting body. 
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1. Introduction 

The tidal dynamics of inspiralling compact binaries (involving neutron stars and/or 
black holes) has been the subject of vigourous investigation in the last several years, 
motived by the exciting prospect of measuring tidal signatures in the gravitational 
waves emitted by such systems. Some of this work has focused on calculating the 
influence of the tidal coupling on the gravitational waves, and estimating the accuracy 
with which the tidal deformation of each body can be measured [T1I2|S1I11|S]. Some has 
focused on calculating the tidal deformation of neutron stars in the post-Newtonian 
approximation to general relativity [H [7l [8] and in the full theory [H [9l [10] . And some 
has focused on the tidal deformation of nonrotating black holes [TTJ [T^l [131 [HI \T5\ . 

An issue that is central to all these investigations is the dependence of adopted 
measures of tidal deformation on the coordinates employed to describe the spacetime 
geometry. In the case of neutron stars, the coordinate independence of the relativistic 
Love numbers which measure the tidal deformation of the body's external gravitational 
field was firmly established by Damour and Nagar |S] and Binnington and Poisson 
[TO] . In the case of nonrotating black holes, however, these gauge-invariant Love 
numbers were shown to vanish [TU] , and the identification of nonvanishing, coordinate- 
independent measures of tidal deformation has remained an open problem. For 
example, Poisson and Vlasov [15] rely on light-cone coordinates to describe the 
geometry of a deformed black hole, while Damour and Lecian [2] rely on Weyl 
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coordinates in a context of stationary and axisymmetric tides. Our main objective 
with this paper is to remedy this situation by providing a complete description of 
the intrinsic and extrinsic geometries of a tidally deformed event horizon, and fully 
clarifying the coordinate dependence of all horizon quantities. In particular, we 
introduce meaningful and practical measures of the tidal deformation of an event 
horizon. 

The central assumptions in our work are that the unperturbed black hole 
is nonrotating and described by the Schwarzschild solution to the Einstein field 
equations, and that the tidal deformation is sufficiently small that it can be described 
accurately to first order in a perturbative treatment. Otherwise our formulation is 
completely general: the tide can be either static, slowly varying, or fully dynamical, 
and there is no requirement that it be axisymmetric. Our description of a tidally 
deformed event horizon relies on two major theoretical foundations. The first is a 
Gauss-Codazzi theory of null hypersurfaces embedded in spacetime, an extension of 
the standard theory of (spacelike and timelike) hypersurfaces formulated in terms of 
first and second fundamental forms. This material is developed here ah initio, in 
spite of the fact that similar formalisms are extant in the literature (for example, in 
Refs. [inillTlITH]); our version is presented in a form directly suited to our application 
to perturbed event horizons. The second foundation is a covariant and gauge-invariant 
formulation of black-hole perturbation theory, as summarized in the work of Martel 
and Poisson [T^. 

Our description of a null hypersurface embedded in spacetime is tied to its 
generators, the congruence of null geodesies that trace the hypersurface. We label 
each generator with two comoving coordinates = (a, /3) (with the index A 
running over the values 2 and 3), and we let A be a running parameter on each 
generator. The hypersurface is charted with the intrinsic coordinates (A, a"^), and its 
(degenerate) intrinsic geometry is fully characterized by the explicitly two-dimensional 
metric 7_as, the analogue of the first fundamental form of a (spacelike or timelike) 
hypersurface. The extrinsic geometry, on the other hand, is characterized by a scalar 
K (a generalization of the black hole's surface gravity), a vector uja, and a tensor 
Kab ; these are analogous to the second fundamental form of a (spacelike or timelike) 
hypersurface. We examine how these quantities transform under reparameterizations 
A — >■ A(A, a^) of the generators, and show that while the extrinsic geometry of the null 
hypersurface depends on the parameterization, the intrinsic geometry is independent 
of the parameterization. When applied to an event horizon, this observation becomes 
one of the central results of this paper: the intrinsic geometry of a black-hole horizon is 
invariant under a reparameterization of the horizon's null generators. This statement 
implies that any measure of tidal deformation that derives from the induced metric 
7ab is necessarily invariant under reparameterizations. 

This result can be restated in terms of gauge transformations, small deformations 
x" — >■ x^+f^ of the coordinates employed in the unperturbed spacetime. With regards 
to transformations of the spacetime coordinates x°', the horizon quantities 7ab, wa, 
and Kab are a collection of scalar fields expressed entirely in terms of the hypersurface 
coordinates (A, a"*). As such they are independent of the spacetime coordinates, and 
therefore immune to gauge transformations. As a matter of principle, therefore, all 
horizon quantities are gauge-invariant quantities. The situation, however, is made 
more subtle by a matter of practice, our identification of the generator parameter 
A with the advanced-time coordinate v of the underlying spacetime. With this 
identification, a transformation of the spacetime coordinates is necessarily associated 
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with a reparameterization of the null generators, and the horizon quantities acquire 
a gauge dependence that is inherited from their dependence on reparameterizations. 
In this context, the results summarized in the preceding paragraph can be stated as 
follows: While the extrinsic geometry of a perturbed event horizon depends on the 
choice of gauge, the intrinsic geometry is gauge-invariant. 

Our Gauss-Codazzi theory of null hypersurfaces is developed in Sec. [2j In Sec. [3] 
we examine a nonrotating black hole deformed by an arbitrary distribution of matter, 
describe its geometry in terms of a perturbed Schwarzschild metric, and compute 
the horizon quantities 7a_b, ujai and ICab to first order in perturbation theory. In 
Sec. [4] we specialize our results to tidal deformations produced by a remote distribution 
of matter. Adopting a specific choice of gauge (the "Killing gauge"), we involve 
the vacuum field equations near the horizon to express the horizon quantities in 
terms of the well-known master functions ^'even (the Zerilli-Moncrief function) and 
(the Cunningham-Price-Moncrief function) of black-hole perturbation theory. 
In Sec. [5] we consider two applications of our formalism, the first involving a slowly- 
varying, quadrupolar tidal field imposed on the black hole, the second involving a close 
parabolic encounter between the black hole and a small orbiting body. An appendix 
( Appendix A I contains mathematical developments regarding the late-time behaviour 
of the horizon quantities. 



2. Differential geometry of null hypersurfaces 

To guide the development of a theory of perturbed event horizons it is helpful to 
formulate a differential geometry of embedded null hypersurfaces. The main goal is to 
arrive at a set of Gauss-Codazzi equations that apply to the null case instead of being 
restricted to usual cases of timelike or spacelike hypersurfaces. The developments of 
this section rely on material presented in Sees. 3.1 and 3.11 of Ref. [20j. 



2.1. Generators, vector basis, and intrinsic coordinates 

A null hypersurface is generated by a congruence of null geodesies that are described by 
the parametric equations x°' = a;"(A, a^), in which A is a running parameter on each 
generator, and = (a, /3) are generator labels that stay constant on each generator; 
uppercase latin indices such as A run from 2 to 3. The null vector field 



(2.2) 



is tangent to the congruence of null generators, and 

are spacelike displacements vectors that point from one generator to another. These 
are orthogonal to fc", fc^e^ = 0, and their mutual inner products are 

lAB QapeAe-B- (2-3) 
The definitions imply that the vectors satisfy the Lie-transport equations 

fc'^^e^ = e%pk^ e%pe% = eg^^e^, (2.4) 

in which a semicolon indicates covariant differentiation in spacetime, with a connection 
compatible with ^q,^. The basis is completed with a second null vector TV" that cuts 
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across the hypersurface; its normalization is chosen so that N^k" = — 1, and the vector 
is also required to satisfy iV^eJ^ = 0. 

We select (A, a"*) as intrinsic coordinates on the hypersurface. In the spacetime 
coordinates a displacement within the hypersurface is described by dx°' = k" dX + 
da'^, and the intrinsic line element is 

ds'^ =jABda^da^. (2.5) 

This shows that 7^13 (A, a'^) acts as a metric on the hypersurface. In this generator- 
adapted coordinate system, the degenerate metric is explicitly two-dimensional. We let 
7^^^ denote the matrix inverse to 7ab, and we let be the connection compatible 
with the two-dimensional metric; the associated covariant-derivative operator is 
denoted Va- We use 7^13 and its inverse to lower and raise uppercase latin indices. 



2.2. Gauss-Weingarten equations 

The tangent vector fields admit the following set of Gauss-Weingarten equations: 



k . i^k^ — t\i k ^ 



a P 



= BabN" 



^A ^B 
+ ICABk'' 



-A-pi^ 



^ AB^C 



-B;I3^A- 



These equations define k, uja, Bab, I^ab, and F^^. Explicitly, 



^N^k^.pk^, 



UJA 

Bab 
ICab 



I3^A' 



= ka-pe'Xe'^ = Bba, 



P^B 



ICba, 



(2.6) 
(2.7) 
(2.8) 

(2.9) 
(2.10) 
(2.11) 
(2.12) 
(2.13) 



FcAB = ecae^.^e^ — Tcba, 

where Tcab ■— Jcd^ab- These equations reveal that while jab, Bab, and F^^ 
characterize the intrinsic geometry of the hypersurface, k, uja, and ICab characterize 
its extrinsic geometry. 

Each equation in the set (2.6)-(2.8) can be viewed as an expansion of a vector 
field (defined on the left-hand-side) in terms of the hypersurface basis N"-, fc", and 
6^4. The absence of terms proportional to N" in Eqs. (2.6) and (2.7) is a consequence 
of the fact that fc" is nul l eve rywhere on the hypersurface. The absence of a term 
proportional to in Eq. (2.6 I follows from the identity eAak^pk^^ = — fco-e^.^fc^ = 0, 
which itself follows from the orthogonality of k" and e^. Equality of Bab as defined 

is confirmed by a similar calculation. 



by Eq. (2.7) and Bab as defined by Eq. (2.. 



Equation (2.6 1 states that k" is a geodetic vector field, and k measures the failure of 
A to be an affine parameter. 



The definition of Eq. (2.3 1 and the Gauss-Weingarten equations imply that 
d\^AB 



2B 



AB- 



(2.14) 



It is customary to decompose Bab into irreducible components, 



B 



1 



AB 



^^lAB 



GAB, (2.15) 

with the trace term 8 Bab representing the rate of expansion of the congruence 

of null generators, and the tracefree term gab '■— Bab — 5 ©Tab representing the rate 
of shear. A similar decomposition could also be introduced for ICab ■ 
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The Gauss- Weingarten equations also imply that 

N'^.pk^ = -kN"" + w-^e^, (2.16) 

iV^^e^^ - ^loaN'' + IC/e%. (2.17) 
These equations govern the behaviour of the transverse vector on the hypersurface. 

2.3. Gauss-Codazzi equations 

It is straightforward, following the methods described in Sec. 3.5 of Ref. ^21j, to derive 



from Eqs. (2.6|-(2.8| a set of Gauss-Codazzi equations which express projections of 
the spacetime Riemann tensor in terms of geometric quantities defined on the null 
hypersurface. We have 

R^.^Xc^k^N'^k^el = dxiOA - Oak + B^ujb, (2.18) 
R^,,o.pk^'N''e1e% = ^ a^b - Vb^a - B^JCcb + B^JCca, (2.19) 
R^c..pk''e"AN''e^s = -BxICab - kICab + ^aoob + ^a^b + K.SBcb , (2.20) 

i?Mc^'/3fc''eAfc"4 =-dxBAB + KBAB + BfBcB, (2.21) 

i?Mo/37fc''e^4ec = ^cBab-^bBac-^cBab+ujbBac, (2.22) 
Rfia/ifN'^eAeQel. = ^c^ab - ^ bK-ac + ^c^ab - ^^bICac, (2.23) 

Ra/i-tseAeBeceo = -;^T^{iac1bd ~ IadIbc) + Bac^bd - BadK.bc 

+ ICacBbd-ICadBbc, (2.24) 
where TZ is the Ricci scalar associated with the two-dimensional metric "/ab ■ To arrive 
at these equations we used the fact that the Riemann tensor on a two-dimensional 
metric space can always be expressed as TZabcd = \TZ[iac1bd ~ IadIbc)- We 
also relied on the identity ^codx^AB = ^ aBbc + ^ bBac ~ ^ cBab, which can be 
derived on the basis of Eq. (2.14). 

Insertion of the Gauss-Codazzi equations within the identity 

R^.. = g"^i?aM73. = {-k^Nf - N^k^ + 7^''eAe^)^aM0. (2.25) 
produces the following components of the Ricci tensor: 

Ri^^k^k" = -dx& + Ke- BabB^^, (2.26) 
R^^ak'^e'X = dxuJA - Oak - OaQ + ^bB^ + Qua, (2.27) 
Rajse'XeB = ^[dxK.AB + kICab) - (Va^^b + Vb^a) - 2ujaujb 

- 2{JCfBcB + IC^Bca) + QICab + ICBab + ^Tl-fAB, (2.28) 
where K. := ^^^JCab- 



By involving the Einstein field equations, Eq. (2.261 can be turned into 
Raychaudhuri's equation, 

dxQ = ne- - (TABcr^^ - 8nTafik°'k'^, (2.29) 

with Tapk^k^ representing the flux of matter across the null hypersurface. And by 
extracting the tracefree piece of Eq. (2.21) we obtain an analogous equation for the 
shear tensor, 

dxa^B = - e)f^ B - C^b, (2.30) 
where Cab '■= C^avpk^e.'Ak'' are components of the Weyl tensor. 
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2.4- Reparameterizations 

The geometric quantities ^ab, Bab, uja, and Kab all refer to a selected 
parameterization (A, a'^) of the null generators. We first examine how these quantities 
change under a reparameterization of the form 

A-^A(A,a^), (2.31) 

which represents an independent change of parameter on each generator. The 
differential form of the transformation is expressed as 

d~X^e-^{d\-CAda^), (2.32) 

with 

These are functions of (A, a^) on the hypersurface, and the notation was chosen so as 
to simphfy our expressions below. The inverse transformation is dX = e^dX + ca da^. 

As we saw previously, a displacement on the hypersurface is described by 
dx°' — k" dX + e'^da^, but the reparameterization brings this to the new form 
dx"' ^k°'dX + e'X da^, with 

fc" = e^/fc", = + c^fc". (2.34) 

These vectors have the same interpretation as the old vectors: fc" is still tangent 
to the congruence of null generators, but is renormalized so as to reflect the new 
parameterization, and still points from generator to generator. It is easy to show 
that the new transverse vector must be given by 

iV" ^ e-f' (TV" + ic^c^fc" + c^e^) (2.35) 

to satisfy its defining relations. The inverse transformations are fc" = e~'^fc". 



= - c-'^caA:", and iV" = e^A^" + ie-^c^c^fc" - ^ ^a 

The reparameterization produces the following changes in the geometric 
quantities: 

lAB = lAB, (2.36) 

Bab - el'BAB, (2.37) 
K = e'3(K + a^/3), (2.38) 
UJA ^^A- B/cb + KCA + caOxP + dAP, (2.39) 

K-AB = e"'^(/CAB + i^ACB + ^bca + KCaCb + CbOxCa + Vbca 

+ \c^ccBab - B^ccCB - B^ccca) ■ (2.40) 

The term CBdxCA + ^ bca in the last equation is not manifestly symmetric in the pair 
of indices AB. With the definitions of Eqs. (2.33), however, we find that this can be 
expressed in the form 

„ S^A d^X dX 2fl dX 



riad'^X dX dX 

which reveals the required symmetry. An additional change produced by the 
reparameterization is f = + B^cb + B^ca — Babc^ ■ 
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In the case of infinitesimal transformations described by A = A + (5A(A, a^), 
the partial derivatives are captured by (5/3 :— —d\SX and Sca '■— —OaSX, and the 



transformations of Eqs. ( 2.36 1-( 2.40) simplify. For the purposes of an application 

of the formalism presented below, we assume that the geometric quantities can be 
expressed as 

lAB =1ab+Sjab, (2.42) 

Bab^SBab, (2.43) 

K = Ko + Sk, (2.44) 

UJA = SujA, (2.45) 

ICab = ICab + ^K^AB, (2.46) 

where the "background quantities" kq, and /C^iB assumed to be A- 

independent, and where S^ab, SBab, Sk, Sloa, and SICab are A-dependent 
"perturbations." In this restricted context the transformations reduce to 

SlAB =S-fAB, (2.47) 

6Bab = SBab, (2.48) 

Sk =Sk + dxSl3 + kqS^, (2.49) 

SuiA = Su!A + QaSP + kqSca, (2.50) 

SICab = SICab-IC%S(3 + VbSca. (2.51) 

In the last equation the covariant derivative Vb is evaluated with a connection 
compatible with the background metric Jab- 

We next examine the possibility of transforming the generator labels. A general 
transformation of the form — >■ a"^(A, a^) is excluded, because the dependence upon 
A would imply that is not constant on each generator, in violation of its defining 
property. The remaining freedom is a rigid transformation of the form — ^ a^{a^), 
upon which scalars such as k remain invariant, while tensors such as uja and ^ab 
transform in the usual way. In particular, for infinitesimal transformations of the 
form = + Sa^ , the metric tensor transforms as 

lABipF) = -iab{cP) - ^aScxb - VsfoA, (2.52) 

where Va refers to 7ab, and Sola = ^abScx^ . In this formulation the original metric 
is expressed as a function of the new coordinates (instead of the original coordinates) , 
and the transformation takes the standard appearance of a gauge transformation. 



3. Deformed black hole 



We consider a nonrotating black hole perturbed by a distribution of matter. The 
perturbation is sufficiently small that we can describe it within linearized perturbation 
theory, and to achieve this we rely on the formulation of the theory provided in 
Ref. [IH]. The matter is either flowing across the event horizon, in which case the 
perturbation is sourced by matter, or it is situated outside the black hole's immediate 
neighborhood, in which case the perturbation is in vacuum and describes a tidal 
deformation of the black hole. 
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3.1. Spacetime metric 

The metric of the unperturbed spacetime is Schwarzschild's solution expressed in 
Eddington-Finkelstein coordinates, 

glp dx^dx^ = -/ dv'^ + 2 dvdr + {dO'^ + sin^ 9 dcj)'^) , (3.1) 

with / := 1 - 2M/r. We let x°- = {v, r) and 6^ = {9, 4>). The metric on the unit two- 
sphere is VLab d9^d9^ — d9^ + siv? BdcjP', and its inverse is denoted Sl"^^; covariant 
differentiation compatible with VLab is denoted Da. 

The metric perturbation is denoted Pap, and it is decomposed in tensorial 
spherical harmonics (as defined in Ref. [19]). In the even-parity sector we have 



Pab ^ hab{v,r)Yi9^), (3.2) 

PaB = ja{v,r)YB{9^), (3.3) 

PAB = r^K{v,r)nABY{9^)+r^G{v,r)YAB{0^), (3.4) 

with Y{9^) denoting standard spherical-harmonic functions, Ya DaY, and Yab '■= 
[DaDb + ^({i + l)ilAB]Y- In the odd-parity sector we have 

Pab = 0, (3.5) 

PaB = ha{v,r)XB{9^), (3.6) 

PAB ^ h2{v,r)XAB{0^), (3.7) 

where Xa ■— —e^DbY and Xab ■— \{DaXb + DbXa), with eab denoting the 



Levi-Civita tensor (with component sg^ = sin 9) on the unit two-sphere. The tensorial 
harmonics Yab and Xab are both symmetric and tracefree. The spherical-harmonic 
labels £m are suppressed, and so is summation over these labels. The complete metric 
of the perturbed spacetime is gap — g^p + Pa/s ■ 

Under an even-parity gauge transformation generated by the vector field fa = 
r]a{v,r)Y and fA = r'^ri°'^°'^{v,r)YA, the perturbation fields change according to 



Ahyy = -2dvr]y + -^rjy + —^Vr, (3.8) 

Ahyr = -drTjy - dy7]r " ^f?r, (3.9) 

Ahrr — —2drr]r, (3.10) 

Aj, = -r^S.r;''™'^ - ?7„, (3.11) 

Aj, ^ -r^drv'^^^ - Vr, (3.12) 

AK =-^r,r--vv+e{e+W'', (3.13) 

r r 

AG =~2?7<=™". (3.14) 

Under an odd-parity gauge transformation generated by the vector field fa=0 and 
fA = r^r]°'^'^{v,r)XA, the perturbation fields change according to 

Ahy = -r^dyrj°'^'^, (3.15) 

AK = -r^drrf'^'^, (3.16) 

A/i2 = -2r^'n°'^'^. (3.17) 

These transformations will play a role in the forthcoming developments. 



Intrinsic and extrinsic geometries of a tidally deformed black hole 



9 



3.2. Deformed horizon 

The description of tlic deformed horizon rehes on the geometrical methods reviewed 
in Sec. [2j The event horizon is traced by its null generators, which are identified by 
constant labels a"^ = (o^,/?); we use A = w as a running parameter on each generator, 
and (?;,a^) forms a system of intrinsic coordinates on the horizon. The parametric 
equations that describe the horizon's position in the unperturbed spacetime are v — v, 
r = 2M, and 9^ = a^. In the perturbed spacetime we have instead 

v = v, r = 2M[l + B{v,a'^)], 6*^ = + a^), (3.18) 

where 2MB and S'^ are the components of a Lagrangian displacement vector. This 
vector takes the horizon point identified by (w, a^) in the original spacetime to a point 
also identified by (w,a^) in the perturbed spacetime. We express the displacement 
fields as 

B = biv)Yia^), E^^n^^[r'''\v)YB{a^)+e''''{v)XBia^)], (3.19) 

in which i}^^ is expressed in terms of the intrinsic coordinates . As previously we 
suppress the £m labels, as well as summation over these labels. 

The parametric equations (3.18) allow us to calculate the basis vectors 

In this section we place brackets around a basis index (which refers to the intrinsic 
coordinates a^) to distinguish it from a coordinate index (which refers to the spacetime 
coordinates 9"^). Explicitly, 

F = 1, (3.21) 

k'' =2MdyB ^2MbY, (3.22) 

fc-^ = a^S-^ = [r™">s + C'^'^Xb] , (3.23) 

in which an overdot indicates differentiation with respect to v, and 

e(A)=0, elA)^2MdAB, ef^) = + a^S^. (3.24) 

The null condition /ccfc" =0 gives rise to the first horizon equation, 

b- 4:Mb = hyy{v,2M), (3.25) 

and the conditions fc^e"^^ = give rise to a second set of horizon equations, 

= -(2M)-2 [j,{v, 2M) + 2Mb{v)\ , (3.26) 
^odd ^^(2M)-^K{v,2M). (3.27) 

These equations, along with appropriate choices of boundary conditions, fully 
determine the description of the deformed horizon. 

The basis can be completed with a transverse vector N'^ that satisfies the relations 
NaN" = 0, Nak" = —1, and TV^e"^^ =0. A simple computation reveals that the 
components of this vector are given by 

iV" = ^KrY, (3.28) 
iV =-l + hyrY, (3.29) 
N"^ = {2M)-^n^^ {jrYB + KXb) , (3.30) 
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where all perturbation fields are evaluated at r — 2M. The covariant components of 
the vector are = —1, N,, — ~}^h„Y , and Na — 0. 

To identify the correct solutions to the horizon equations we imagine first an 
artificial situation in which the perturbation is switched off at times larger than vi. 
The spacetime for u > t^i is described by the Schwarzschild metric, and for these 
times the event horizon is correctly identified with the hypersurface r — 2M. To 
locate the event horizon at times v < vi we must smoothly extend r — 2M backwards 
in time, to a null hypersurface in the perturbed spacetime. This surface is described 
by Eq. (3.18), with b{v) restricted to vanish for v > vi. The appropriate solution to 



Eq. (3.25) is therefore 



,-ko(u'-d) 



hyy{v',2M) dv', 



(3.31) 



where kq (4M)^^ is the surface gravity of the unperturbed black hole. The upper 
limit of integration was extended to w = cxd because, by the stated assumptions on 
the perturbation, h^^ is zero in the interval vi < v' < oo. At this stage, however. 



the artifice can be removed and Eq. (3.31) be adopted as the appropriate solution to 
Eq. (3.25) even when the perturbation does not switch off at v — vi. The perturbation 



must still fall off sufficiently fast that the integral converges, and under these conditions 
b{v) will approach zero as u — oo. Because Eq. (3.31) reflects a choice of final 



condition, it is known as a ideological solution to the horizon equation. 



The teleological solutions to Eqs. (3.261 and (3.27) are 



[jy{v' ,2M) + 2Mb{v')\ dv', 
K{v\2M)dv'. 



(3.32) 



(3.33) 



The behaviour of the horizon generators in the perturbed spacetime is now completely 
determined. The solutions to the horizon equations imply that in general, the event 
horizon leads the perturbation by a time interval of order Kg ^ = 4M. 



3.3. Horizon's intrinsic geometry 

As described in Sec. [2] the intrinsic geometry of the event horizon is characterized 
by the induced metric 7ab, which is expressed in the intrinsic coordinates [v^a-^) 
attached to the null generators. According to Eq. (2.14), the w-derivative of the 
induced metric satisfies 

dvlAB = '2Bab = QlAB + 2crylB, (3.34) 

and this equation defines the expansion scalar Q and shear tensor ctab associated with 
the congruence of null generators. The expansion, in particular, can be computed as 



e 



'^9u7, where 7 :— detl'-fAB] 



A computation of the horizon metric involves the substitution of Eq. (3.24) into 



Eq. (2.3). The computation must account for the fact that while the spacetime 
metric is expressed in terms of the coordinates [v^r^O^), the horizon metric will be 
expressed in terms of the intrinsic coordinates (f,a"^). A piece of the computation 
that requires some care involves f7^B(0"^), which must be written as $1^5(0;"^ + S"^) = 
Q.AB + dc^AB^ with the right-hand side expressed in terms of a^. With this 
accounted for, we find that the horizon metric is 

lAB 



{2MfnAB + {2Mf{2BaAB + ^bcDa'EP + ^acDb^^) +Pab- 



(3.35) 
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= {2Mf{^AB+l 

where 



^tracco V 



Yab + I^'^'Xa 
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(3.36) 



7 



3dd 



= 2b{v) - £{i + l)r™"(w) + K{v, 2M), 
= 2C"™"(f) + G(w,2M), 



(3.37) 
(3.38) 
(3.39) 

The square root of the metric determinant is given by y/j = (2Af)^ sinQ!(l _(_ ^ytfaccy-j^ 
with a denoting the intrinsic polar angle on the horizon. 

It follows from these equations that the expansion scalar is 

tracc^ 



1 



9 = 7 

while the shear tensor is 

The expressions for j^'^^^'^^ 
Eqs. ( |3^ -( [3^ . We obtain 



(2M)^ 7°™"r^B + 7 



odd 



Xai 



(3.40) 



(3.41) 



and 7° can be simplified with the help of 



odd 



7^ 



(2M)-i{[^(£ + 1) + l]biv) - K,{v,2M)+t{t+l){2M)-^j,{v,2M) 
+ {2M)dyK{v,2M)Y 

(2A/)-i|-26(w) - 2{2M)-^j^{v,2M) + 2Md^G{v,2M)Y 

{2M)-^i^-2{2M)-^h^{v,2M) + (2M)"i9„/i2(w, 2M)|. 

The Ricci curvature associated with the metric of Eq. (3.36) is given by 
1 



2M2 



1 



1 



1+ o(^-l)(^ + 2) 7*''^'='= + -^(£+1)7"=™" Yia"^) 



(3.42) 
(3.43) 
(3.44) 

(3.45) 



This indicates that the metric's geometrical information is contained within + 
i£(£ + 1)7°™'! = 2b + K + + l)G; the remaining information is entirely about 
the choice of intrinsic coordinates, in particular, the fact that they are attached to 
the horizon's null generators. To flesh out this last point we recall that, according to 
the discussion near the end of Sec. 2.4 the freedom to change the generator labels a"^ 
is limited to a u-independent rotation of the form — > a'^{a ^). F or infinitesimal 
changes — + 6a^ the transformation is described by Eq. (2.52 1. If we choose 

da^ = r!^^ {C'Yb + C'^'^Xb) , (3.46) 



where (^°'^™ and C° are constants, we find that 7 



7^ 



and 7°'^'^ change 



according to 

^traco ^ ^traco _^ ^ 1)(<=™", (3.47) 

^cvcn ^ ^cvon _ 2<"=™", (3.48) 

^odd ^ ^odd _ 2C<=™". (3.49) 

We observe that the combination + + i)'y<=™n is unaffected by the 

transformation, which confirms its role as carrier of geometric information. At any 
given time (but at only one such time), and C"^*^ can be chosen so as as to make 
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^ovcn ^odd yanisli. At this time, say v = vq, have that the horizon metric is 
given by 

7AsK,«^) = (2M)2r!^B[l + 7'— (3.50) 



with 7 



trace , _,tracc 

■ /new 



,cvcn 
Id ' 



7'^^^'=(i;o) = 2b{vo) + K{vo, 2M) + -£{£ + l)G(«o, 2M). 



(3.51) 



At other times v ^ vq, the tracefree terms proportional to Yab and Xab will no longer 
vanish, and the metric will return to its general form of Eq. (3.36). 



3.4- Horizon's extrinsic geometry 

The horizon's extrinsic geometry is characterized by k, uja, and /Cab, as defined by 
Eqs. (2.6|-(2.8). Computation reveals that 

K = Ko(l + kY), 



where kq = (4M) ^ is the unperturbed surface gravity, and 
k = - {2Mdrh^v - AMdyKr + K,. + 26). 



We also get 
with 



'YA+iO°'"'XA 



Ddd 



^hyr - ^drjv + (2Af) '^jy + ^dyjr + b, 



And finally, we get 
Kai 



^2AmAB + IC^'^'^riABY + K^^'^'^Yab 



jL^odd V 



with 



^trace 
^odd 



{2MfdrK - 2MK - 2Mb, 



e{l+l)2M£, 
1 

" 2 
_4Af ^°dd ^ 



2Af/i„ 



1 



^(^ + 1)> 



{2MfdrG - 2MG, 

drh2- 



(3.52) 
(3.53) 
(3.54) 

(3.55) 
(3.56) 

(3.57) 

(3.58) 
(3.59) 
(3.60) 



In these expressions, all perturbation fields and their derivatives are evaluated at 
r = 2M, and all spherical harmonics are expressed as functions of a^. We note that 
the computation of /Cab requires the same level of care as the previous computation 
of 7ab: the unperturbed expression — j'Oab must be evaluated at r = 2Af (1 + B) and 
6^ — + S"^ and combined with the terms that arise from the metric perturbation. 
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3.5. Gauge transformations 

We next work out how the various horizon quantities introduced previously are affected 
by a gauge transformation of the form 

x" ^x^ + f, r = r^^iv, r)Y{e^) (3.61) 

and 

e^^e^ + f^, f^ = VL''''[r,''^-^{v,r)YB{e^) + il°''\v,r)XB{e^)\. (3.62) 

We also have that /a = r^^Asf^ = (r2?7'=™")YA + {r^r]°'^'^)XA. The gauge 
transformation affects the coordinate description of the horizon. Recalling Eq. (3.18), 
we find that 5, and i^"'^'^ change according to 

A6 = (2M)-^/„(w,2M), (3.63) 
Af™" = ?7'=™"(w,2M), (3.64) 
A^°dd ^^°dd(^ 2M). (3.65) 

A complete listing of the corresponding changes in the metric perturbation can be 
found in Sec. 13.11 

With these rules it is easy to show that the quantities associated with the horizon's 
intrinsic geometry change according to 

These results imply that 7_ab, ©, and uab are all gauge invariant, and we conclude 
that the horizon's intrinsic geometry is gauge invariant. This is not a surprising 
conclusion. The intrinsic metric is a collection of scalar fields with regards to 
transformations of the spacetime coordinates x°' , and it is expressed entirely in terms 
of the intrinsic coordinates (A, a^). As such it is as a matter of principle immune to 
a gauge transformation. The fact that A is identified with the spacetime coordinate 
V adds a small complication to this argument, because ^ab could in principle be 
sensitive to a change in v. The identification associates a gauge transformation 
on z; to a reparameterization of the generators, as was described in Sec. |2.4[ But 

= (2M)^f2yiB, the induced metric on the unperturbed horizon, is independent 
of V, and the results displayed in Eqs. ( 2.47 1-( 2.51 1 reveal that an infinitesimal 
reparameterization has no effect on 5^AB^ the metric perturbation. The conclusion, 
therefore, remains valid regardless of the identification X = v. As an additional remark, 
we recall that the invariance of ^^ab under general (large) reparameterizations was 
established in Eqs. ( |236l )-( |2lo| . 

On the other hand, the quantities associated with the horizon's extrinsic geometry 
change according to 

A/c = -AMdy {dyrjr + Korjr) , (3.67) 

^^cvcn =-(dyr]r + K0Vr), (3.68) 

^^tracc ^ _^2MdyT]r + ^^{i + l)Vr, (3.69) 

A/C°™" = -7]r, (3.70) 

and 

Aa;°'*'^ = = A/C°^^ (3.71) 

In the even-parity sector, the changes in the extrinsic geometry are all associated 
with rir{v,2M) = rj'" {v,2M), which describes a change in v; there are no changes in 
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the odd-parity sector. As before we can observe that since k, uja, and Kab are all 
spacetime scalars expressed entirely in terms of the intrinsic coordinates, they should 
all be immune to a gauge transformation. But as before we can identify a change in 



V with a reparameterization of the generators, and infer from Eqs. (2.47|-(2.51 1 the 



effect of the reparameterization on the ext rinsic g eome try. With 5\ identified with 
?7r(w, 2M)y(a'^), we quickly find that Eqs. (2.47)-(2.51 ) reproduce the statements of 
Eqs. ( |3!67| -( [3JI] ). 



It is easy to identify four linearly-independent quantities, formed from fc, w''™", 
^tiace^ and /C°™", that are invariant under infinitesimal reparameterizations. We 
choose 



■01 := w 
ip2 := K 



04 := /C*''""'= + I [£{£ + 1) + 1] /C°™" - 



2 ' 
2Mcj'=™". 



(3.72) 
(3.73) 

(3.74) 
(3.75) 



The first three combinations can be shown to be pieces of the spacetime Riemann 
tensor evaluated on the deformed horizon. Indeed, inserting the results obtained in 
Sec. [O] within Eq. ( |2.18[ ) yields 



Kok)YA 



•odd 



Xa- 



(3.76) 



It is easy to show that the left-hand side is invariant under infinitesimal 
reparameterizations, and this guarantees that ipi and 0;°'^'^ also must be invariant. 



Similarly, we find from Eq. (2.28) that 



1 



/C*''""° + Ko/C""=" -I- ^£{1 + l)w""™ - ^k 

+ ]£{£ + 1)7' + W + + 2)7' 

4 o 



i^odd odd 



M7°^< 

odd 



^abY 
Yab 



M7' 



-1°'"']Xab: 



(3.77) 



invariance of Rap^A^s ^^'^ Iab under infinitesimal parameterizations guarantees that 
il)2 and 03 also must be invariant. The fourth quantity, ip^, does not appear to be 
related in a similar way to a piece of the spacetime Riemann tensor. 



4. Tidal deformations 



The formalism developed in the preceding section is very general, and it can 
accommodate black-hole deformations created by matter flowing across the event 
horizon, or by matter situated outside the black hole's immediate neighbourhood. 
The formalism is also general relative to the choice of gauge, because the relations 



between the horizon quantities (such as 7* 



7" 



7' 



odd 



UJ 



odd 



/C™", and JC°'^'^) and the metric perturbation are valid in any gauge. How the horizon 
quantities change under gauge transformations (or better stated, reparameterizations 



of the horizon's null generators) was described in Sec. 3.5 
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In this section we specialize the situation to a tidal deformation of a black hole 
created by a remote distribution of matter. We incorporate the vacuum field equations 
into our analysis to relate the horizon quantities to the well-known master functions 
\l/even and 4'odd of black-holc perturbation theory (defined below) . We next introduce 
a geometric notion of tidal displacement on the event horizon, and describe how the 
tidal bulge is related to the applied tidal field. 



4-1. Master functions 

Gauge-invariant definitions of the master functions were provided in Ref. |19) . In the 
even-parity sector, 4' oven is the Zerilli-Moncrief function [221 El] defined by 



A 



where A := i{l +1) = /i + 2, /i := (£ - 1)(^ + 2) = A - 2, and where hab := 
hah - Va£6 - VbEa, K := K + \\G - 2r°-£a/r, with ja - \r'^^aG, are gauge- 

invariant combinations of metric perturbations. We use the notation ;= dr/dx°'^ 
Vq is the covariant-derivative operator compatible with the two-dimensional metric 
g^j^dx'^dx^ = —fdv^ + 2dvdr, and as usual the spherical-harmonic labels £m are 
omitted. The Zerilli-Moncrief function is known to satisfy the Zerilli equation |22) . 
which is a two-dimensional wave equation with an effective potential and a source 
term constructed from the energy-momentum tensor of the matter distribution. 

In the odd-parity sector, is the Cunningham-Price-Moncrief function p4l[25] 
defined by 

*odd := -e"" hah'r - -rah'r) , (4.2) 
/i V r / 

where Sab is the Levi-Civita tensor on the two-dimensional manifold with metric 
and ha '■— ha — \y + i'ah2/r is a gauge- invariant combination of metric 
perturbations. The master function is known to satisfy the Regge- Wheeler equation 
|26| , another two-dimensional wave equation with an effective potential and a source 
term. The Regge- Wheeler equation is also satisfied by another choice of master 
function, the original Regge- Wheeler function [21]; in vacuum this is equal to half 
the time derivative of the Cunningham-Price-Moncrief function. 

4-2. Killing gauge 

To relate the horizon quantities to ^'ovcn and 4'odd it is convenient to adopt a "Killing 
gauge" defined by 

Pcptl^ = 0, (4.3) 

where t"' is the timelike Killing vector of the Schwarzschild spacetime. In the 
coordinates {v, r, 6, 4>) we have that t" = (1, 0, 0, 0), and the gauge conditions translate 
to 

hvv = hyr =^ jv ^0 (4.4) 

in the even-parity sector, and 

hy = (4.5) 

in the odd-parity sector. These conditions apply in a neighbourhood of the event 
horizon. 
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An immediate virtue of the Killing gauge is that it preserves the coordinate 
description of the event horizon, which continues, even in the perturbed spacetime, 
to be described by r = 2M and 9"^ = a"^. In the terminology of Poisson and Vlasov 
|15] . the Killing gauge is a horizon-locking gauge. This can be seen at once from 



Eqs. (3.31 1-( 3.33 1, which imply that 

todd _ 



b{v) = r 







(4.6) 



whenever h^y = jy ^ hy = at r = 2M. We remark that while the light-cone gauge 
adopted by Poisson and Vlasov also has the property of being a horizon-locking gauge, 
the Killing gauge adopted here is quite distinct from the light-cone gauge. 



4-3. Near-horizon analysis 

To calculate the horizon quantities we must integrate the perturbation equations 
in a neighbourhood of the event horizon; these are listed in Sees. IV B and V B 
of Ref. [19]. In the even-parity sector this can be accomplished by inserting the 
expansions hrr — ho{v) -\- hi{v){r — 2M) -\- h2{v){r — 2M)^ + •••,> = jo{v) -\- ji{v){r — 
2M) j2{v)ir - 2Mf + ■■■, K = Koiv) + Ki{v){r - 2M) -f K2{v)lr - 2MY + ■■■, 
and G = Gq{v) + Gi{v){r - 2M) + G2{v){r - 2MY + • • ■ within the perturbation 
equations, and solving order-by-order in r — 2M. Such an analysis reveals that jo(''') 
is unconstrained by the field equations, that Kq(v) = and Gq{v) = 2ii~^djo/dv, and 
that Ki{v) and Gi{v) must satisfy the differential equations 

■ KoKi ^ -4AKgjo (4.7) 



and 



dv 

dGi 
dv 



kqGi = 



4ko 



dv'^ 



djo „ 2 ■ 



(4.8) 



where kq := (4M)~^ is the surface gravity of the unperturbed horizon. 

The fact that the differential operator acting on Ki and Gi is d/dv -\- kq, instead 
oid/dv — Ko as in Eq. (3.25), implies that one should not look for teleological solutions 
to these equations: the presence of e'^°'" instead of e~''°" within the integrals would 
prevent them from converging if the integrations were unbounded. We work instead 
with the most general solutions 



Ki{v) = Ki{vo)e-^'>^''-'">^ - AXh 



~Koiv~v') 



jo{v')dv' 



(4.9) 



and 



Gi{v) = {Giivo) + 



4ko 



dv 



(vo) - (A - l)KQjo{vQ 



4ko 



djo 

dv 



(i;)-(A-l)«;oJo(«)-(A- 




-Ko{v-Vo) 



e-''»(''-''")jo{v')dv' 



(4.10) 



in which the initial values Ki{vo) and Gi{vo) are not determined by the requirements 
that Ki{v — > oo) and Gi{v — oo) — >■ 0. (We explore these issues further in 



Appendix A below.) We do not need expressions for ho{v), hi{v), ji{v) (nor other 
coefficients in the expansions), and the master function can be shown to be given by 
*cvon = ifJ'Ko)~'^djo/dv at r = 2M. 
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In the odd-parity sector we substitute the expansions hr ~ hrQ{v) + hri{v)(r — 
2M) + hr2{v){r-2M)^ + --- and /12 = h2o{v) + h2i{v){r-2M) + h22{v){r~2M)^ + --- 
into the perturbation equations, and solve order-by-order in r — 2M. Such an 
analysis reveals that hro{v) is unconstrained by the field equations, that /i2o(^) = 
{2^KQ)~^dhrQ/dv, and that /i2i(w) satisfies the differential equation 



dv 

The general solution is 
h2i{v) = I h2i{vo) 



d/i2i , , 1 

+ Koh2l = 



AKq — h IlIKgh,. 



dv 



1 



dhj-Q 
dv 



(vq) + (A - l)Ko^ro(wo) 



rO 



-Ko(t)-iio) 



(4.11) 



1 



^{v) + {X-l)Kohroiv) + {X-3)Kl I e-"° /^,oK)rf«' 



(4.12) 



and the master function can be shown to be given by 5'odd = ^(Af^o) ^dhro/dv at 
r = 2M. 



Horizon quantities 

With the results obtained in the preceding subsection, we find that the horizon 
quantities defined in Sees. |3.3| and |3.4| are given by 

,trace 



7 

7° 

i 

k 



odd 



j^tracc 
^cvcn 



odd 



0, 

2ko*cvc„(v,2M), 
-2Ko*odd(t',2M), 

0, 
1 

-//KO*cvcn(w, 2M), 

-^fiKo'fodd{v,2M), 



^(t;o) + -*cvc„(wo,2M) 



-K;o(u-t)o) 



Av,2M) 



+ ^(A - 3)ko / e-"''("~''''*evcn(«', 2M) di;'. 



^ («o)- 2*odd(fo,2M) 



'"°^"-"''^ + ;^*odd(t',2M) 



(4.13) 
(4.14) 
(4.15) 
(4.16) 

(4.17) 
(4.18) 
(4.19) 



(4.20) 



(4.21) 

where /C*''^'=°(uo),/C'=™" (wo), and /C°^'i(uo) can be expressed in terms oiKi{vo), Gi{vo), 
and joM- We recall that X -.^ i{i + 1) = fi + 2 and fi := {i ~ l){e + 2) = A - 2. 

The results display a pleasing symmetry (up to signs, which are inherited from the 
definitions of the master functions) between the even-parity and odd-parity sectors. 
In the case of the intrinsic-geometry quantities and 7°'^'^, the symmetry is gauge 
invariant; in the case of the extrinsic-geometry quantities a;''™" and uj°'^'^, ^cvcn ^^^^ 
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KP'^'^ , the symmetry is a property of the KiUing gauge adopted here (it is not, in 
particular, a property of the hght-cone gauge [H]). Another remarkable property of 
the Killing gauge is the fact that fc = 0, so that the surface gravity of the perturbed 
black hole is k = kq = (4M)~^ 

The expressions for /C'''^'^° and /C''™" given previously were simplified relative to 
the more primitive expressions obtained in terms of jq. These, however, involved the 
combination 

JO- no r e"'^^''''^Joiv')dv' ^ joMe-^"^^-""^ + P e'^^^-^'^p^ dv' , (4.22) 

Jvo Jvo dv' 

which could readily be expressed in terms of ^'cvcn — {lJ'i^o)~^^djt)/ dv. A very similar 
simplification was achieved in the case of JC°'^'^. 

The gauge-invariant quantities defined by Eqs. (3.72)-(3.75l are easily shown to 
be given by 

V'l = ^M'«o<9„*cvcn(w,2M), (4.23) 

4'2 = -^a.*cvo„(w, 2M) - Ko*cvc„(i', 2M), (4.24) 

V'3 = -^A^Ko*cvc„(w,2M), (4.25) 



V'4 



-K.a(v-vo) 



JC'^^'^vo) + ^(A + l)/C'=™"(i;o) + J(A + l)*cven(^'o, 2Af) 
- J(2A- l)^evcn(«,2M)- i(A/i + 3)Ko / e-«°('^-''')^evcn(f',2M)(i«'. (4.26) 



4-5. Intrinsic geometry and tidal displacement 

The results obtained in the preceding subsection imply that the induced metric on the 
event horizon simplifies to 

JAB ={2Mf{^lAB + 5-fAB), (4.27) 

5-iAB = 2ko [*cvc„(w, 2M)Yab - ^odd{v, 2M)Xab] (4.28) 
in the case of a tidally deformed black hole. From Eqs. (3.40) and (3.41 1 we also get 
9 = (4.29) 

and 

aAB = M[9„^'evcn(«, 2M)Yab - 5„*odd(w, 2M)Xab] ■ (4.30) 

The fact that the expansion vanishes to leading order in perturbation theory can 
be derived directly from Raychaudhuri's equation: The reduction of Eq. (2.29) to 
vacuum and to first-order perturbation theory is d^Q ~ kqQ, and this implies (with 
an appropriate choice of final condition) that O must vanish. 

The reduction of Eq. (2.30) gives an expression for the Weyl tensor evaluated on 
the event horizon: 

Cab = (kq - dy)aAB, (4.31) 

where Cab '■= C ^avpk^ e'^^A)^^ ^^(b)- This equation can be integrated to relate the shear 
tensor to the Weyl tensor; the appropriate teleological solution is 

poo 

aAB{v,a^)^ e--'>^-'~-^CAB{v',a^)dv'. (4.32) 
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This equation implies that the shear tensor anticipates the behaviour of the Weyl 
tensor by a time interval of order Kq ^ = 4M. If the Weyl tensor is identified with the 
tidal field acting on the black hole, and if the shear tensor is adopted as a measure 
of tidal deformation, then we have the statement that the tide leads the applied field 
by a time interval of order AM. This observation was already made by Fang and 
Lovelace [13] in a more restricted context (and by Hartle [27] in the case of a rotating 
black hole), and we find here that it holds in all generality as a consequence of the 
teleological nature of the event horizon. We remark that in the case of a Newtonian 
body made up of a viscous fluid, the tide would be lagging instead of leading (when 
the body is nonrotating) , and that the time interval would be proportional to Rv/M, 
with R denoting the body's averaged radius, v its kinematic viscosity, and M its mass. 

Another meaningful measure of tidal deformation comes from the Ricci curvature 
scalar associated with the metric of Eqs. (4.271 and (4.28 1. This is given by Eq. (3.45 1 
with 7"'^'^° 



and 7' 

n = 



1 

2M2 



2ko*ovc„(w,2M): 
^ + + + 2)Ko*cvcn(^', 2M)r(a^) 



(4.33) 



It is helpful to convert this into a dimensionless tidal displacement field p{v,a^) by 
identifying TZ with the curvature of a two-dimensional surface embedded in a flat, 
three-dimensional space. We describe this surface in spherical coordinates {r,a^) by 
the parametric equation r = 2M[1 + p{v, a")] with p = e{v)Y{a^), and demand that 
its curvature be equal to TZ. We thus obtain 2]VPn =[l + {e-l){i + 2)eY{a^)], and 
the identification 

p{v, a^) = ^e{£ + l)Ato«'evcn(«, 2M)Y{a'^) (4.34) 

follows immediately. Once more summation over the omitted spherical-harmonic 
labels im is understood. 

The absence of trace terms in Eqs. (4.281, (4.30), and (4.31) implies that each 
tensor S^ab, o'ab, and Cab possesses only two independent components. Introducing 
the basis vectors 

a^ = [l,0], (3^ = [1,1/ sin a] (4.35) 
{a, /3), we take the independent components of S^ab 



in the horizon coordinates a 
to be 



7+ 



7x 



AB, 



a 



^(3^ + p-^a^)<57^B 



(4.36) 
(4.37) 



The independent components cr+^x and C+.x of the shear and Weyl tensors are defined 
in a similar manner. These quantities are closely analogous to the gravitational-wave 
polarizations /i+.x that can be defined in the wave zone of an asymptotically-flat 
spacetime. 



5. Applications 

5.1. Slowly-varying quadrupolar tidal field 

As an application of the general formalism developed here we revisit the situation 
examined by Poisson and Vlasov [TS], that of a black hole deformed by a slowly- 
varying tidal field. To simplify the discussion we neglect the nonlinear terms included 
in Ref. [H], and we specialize the tidal field to a pure quadrupolar form. 
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As described in Sec. II of Ref. [15], the black hole's tidal environment is described 
by the tidal moments £jk{v) and Bjk{v). These quantities are symmetric-tracefree 
(STF) Cartesian tensors that represent the components of the spacetime Weyl tensor 
evaluated far away from the black hole; latin indices j and k (and so on) run over the 
values 1, 2, and 3. The tidal moments give rise to the tidal potentials 

= EpqnPn", (5.1) 
= P/£pkn\ (5.2) 
= 2P/P,''£p, + P,k£' (5.3) 



and 



b; ■.^e.p.nPB^n^, (5.4) 

B;, e,p,nPB\Pl + ekp,nPB\P-^, (5.5) 

where the label "q" stands for "quadrupolar," il^ :— [sin a cos /3, sin a sin /3, cos a] 
is a Cartesian unit vector constructed from the generator labels — {a, (3), 
Pjk ■= Sjk — flj^k is a projection operator to the subspace transverse to , and ejkn 
is the Cartesian permutation symbol. The vector potentials S'^ and B'^ are transverse, 
in the sense that E^il^ = = B'^ft^ . In addition to being transverse, the tensor 
potentials fj^, and S^^, are also tracefree, in the sense that S^^E^j^ = = S^'^Bj/,. In all 
manipulations involving Cartesian tensors, latin indices are lowered and raised with 
the Euclidean metric 6jk- 

The vectorial and tensorial potentials can be converted to angular components 
by means of the transformation matrix il'^ :— dil,^ /da^. We thus introduce 

£\:^£P'^, £\s:^£],^l\^% (5.6) 

and 

B\:^B]n\, B\^:=B;,n^n%. (5.7) 

As shown in Sec. II of Ref. [T5], these angular potentials can be expressed as expansions 
in spherical harmonics of degree £ = 2. We have 

cq _ \ ^ c \r2,m cq _ 1 \ ^ c 1^2, m cq _ \ ^ c ^'i-ra (r o\ 



and 



= J E ^-^1'"' ^AB = E ■ (5.9) 



The sums are carried out from m = —2 to m = 2, the coefficients £m and Bm are 
related to £jk and Bjk and depend on v only; the spherical harmonics are functions 
of a^. These expansions reveal that £jk{v) gives rise to a perturbation of even parity, 
while Bjk{v) gives rise to a perturbation of odd parity. 

Solutions to the perturbation equations corresponding to a black hole deformed 
by a quadrupolar tidal field were constructed by Poisson and Vlasov [l3]. The 
construction assumes that the tidal moments vary slowly, in the sense that the 
timescale r associated with these variations (denoted TZ in Ref. [15]) is very long 
compared with the black-hole mass. The solutions were provided in the light-cone 
gauge, but it is easy from these results to obtain the gauge-invariant master functions. 
The relations, in fact, are the same as in the Killing gauge adopted in Sec. |4.2| we 
have that 4'even(w,2M) = 2MG{v,2M) and «'odd(w,2M) = -(2M)-i/i2(w, 2M), 
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where G and are obtained in the light-cone gauge. Importing the results of 
Ref. [inj — summarized in their Eqs. (6.10), (6.18) and Table XIV — we find that 
G{v,2M) = -^A'PSrn and h2{v,2M) 



so that 



i(u, 2M) = --M^£„,{v) 1 + 0{M^/t^) 



*odd(i',2M) =-Af3S„(w) 1 + 0{M^/t^) 



(5.10) 
(5.11) 



Notice that these expressions involve £m and Bm only, and not their derivatives with 
respect to v, which would contribute fractional corrections of order M/r and {M/tY'. 
As explained in Sec. VI A of Ref. [15], this property results from the freedom to 
redefine the tidal moments according to £jk £jk + PiMSjk + p2M'^Ejk + • • • and 



B,k ^ Bjk + qiMBjk + q2M^Bjk ■ 



where pi, p2, qi, and q2 are arbitrary numbers. 



It is a simple matter to insert Eqs. (5.101 and (5.111 within Eqs. (4.13 )-(4.21 1 



and to calculate the horizon quantities. Because the tidal moments £jk and Bjk vary 
slowly, the integrations can be carried out as in [Appendix A[ by repeated integration 
by parts. After discarding the transient terms that decay exponentially, we arrive at 



7 



0, 



k 

, .even 
Lu 

^ trace 
^cveii 

^odd 



0, 



M^Br, 



M'^Br, 



M''{£„, - 6M£„r + 24M^£m), 



--M^{B,n - 6MB„^ + 2AM^B„ 



(5.12) 
(5.13) 

(5.14) 
(5.15) 
(5.16) 

(5.17) 
(5.18) 
(5.19) 

(5.20) 



These expressions are valid up to correction terms of fractional order (M/t)^; they 
are given in the Killing gauge introduced in Eq. (4.3). 



These results, together with Eqs. (|5.8|) and (|5.9|), imply that in the Killing gauge, 
2 



Qab-^^M^{£\b 



JAB = {2My 
1 

" 4M' 



B 



AB) 



ICab = -2MVLAB + 16M3(f - 4Af£'' + 16Af2f "jr^Ais 



3 



{^AB + B\b) ~ QM{£\^ + B\^) + 24M2 + s^^) 



(5.21) 
(5.22) 
(5.23) 

(5.24) 



as before these expressions are accurate up to terms involving the third derivative 
of the tidal moments. We showed in Sec. |3.5| that jab is gauge-invariant, while k, 
wa, and K-ab are affected by a reparameterization of the horizon's null generators. 
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Figure 1. Orbit of a parabolic encounter between a small body of mass m and 
a black hole of mass M. The orbit's semi-latus rectum is p = 8.1M and its 
eccentricity is e = 1. The particle begins from rest at infinity, reaches a radial 
turning point at r = 4.05M, and returns to rest at infinity. The orbit is displayed 
in a x-y plane constructed in the usual way from the Schwarzschild coordinates 
r and <p, so that x = r cos ip and y = r sin (/>. The coordinates are rescaled by a 
factor of 2M to make them dimensionless; in these units the unperturbed horizon 
(shown in black) is described by a circle of unit radius. The orbital motion is 
calibrated so that (f) = when r = 4.05M. 



Gauge-invariant combinations of these quantities were identified, and in particular we 
have that 



is invariant under infinitesimal rcparamctcrizations; because it originates from a;°'^°" 
and this expression is accurate up to the fourth derivative of the tidal moments. 

5.2. Parabolic encounter 

As a second application of the formalism we consider a parabolic encounter between 
a particle of mass m and a black hole of mass M . We take m to be much smaller 
than M, and we take the motion of the particle to be a geodesic in the Schwarzschild 
spacetime. We give the orbit a semi-latus rectum p — 8.1-M and an eccentricity 
e — 1] the parameterization is such that the radial turning points are situated at 
'"min = + e) = 4.05Af and r„iax = — e) = oo. The orbit has a KiUing 
energy E = m and a Killing angular momentum L ~ 4.0003mM. The particle begins 
from rest at infinity, moves inward, circles approximately twice around the black hole, 
moves outward, and returns to rest at infinity; the shape of the orbit is displayed in 
Fig. [T] Because the turning point is so close to the black hole, the motion is highly 
relativistic when the particle revolves around the black hole, and the tidal interaction 
is highly dynamical. 

We calculate the gravitational perturbations created by the orbiting particle by 
integrating the Zerilli and Regge- Wheeler equations for the master functions ^'cvcn 
and ^'odd- This must be accomplished numerically, and we rely on the time-domain, 
finite-difference code written by Karl Martel; the details of the code are described in 
Refs. [51]. Martel's original code had to be modified to account for a different 




(5.25) 
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Figure 2. Polarization 7+ associated with the intrinsic geometry of a black- 
hole horizon perturbed by a parabolic encounter, calculated at azimuthal position 
/3 = on the horizon's equatorial plane q = ^, which coincides with the orbital 
plane. The polarization is displayed as a function of v/(2M) and is rescaled by 
a factor of m/(2M). All relevant multipoles up to £ = 4 are included in the 
computation. 



choice of odd-parity master function: While Martel's code integrates the Regge- 
Wheeler equation for the original Regge- Wheeler function (which is equal to |c)(5'cvon), 
our modified version of the code calculates instead the Cunningham-Price-Moncrief 
function ^'cvon- The code returns the master functions evaluated as functions of v at 
a fixed radial position r = 2Af(l -I- e) close to the event horizon; in our runs we chose 
e ~ 10-^ 

In Fig. [2] we plot the polarization 7+ associated with the horizon's intrinsic 
geometry, as defined by Eq. (4.36); this is shown as a function of advanced-time v 
at azimuthal position 13 ~ Q on the orbital plane a = ^ ; for this orientation we have 
that 7x = 0. The calculation involves a summation over all multipoles up to (and 
including) ^ = 4; multipoles with £ > 5 give contributions that are too small to be 
visible in the plot. Most of the signal is produced when the particle revolves around the 
black hole, and the plot reveals the rich harmonic structure that a parabolic encounter 
imprints on the tidal deformation of an event horizon. 

In Fig. [3] we plot the polarizations (t+ = ^dv-f+ and C+ — {kq ^ c't,)o-+ of the 
shear and Weyl tensors, respectively; these also are displayed as functions of v at 
position a = ^ and /? = on the event horizon. The figure reveals very clearly that 
the horizon tide (as measured by the shear tensor) leads the tidal field (as measured 
by the Weyl tensor) by a time interval of order = 4M; this feature of the tidal 
dynamics of a nonrotating black hole was discussed in Sec. |4.5| 

Finally, in Fig.|4]we display the shape of the tidal bulge at a selected value of v in 
relation to the position of the orbiting body. The tidal bulge is d escrib ed geometrically 
in terms of the tidal displacement field p{v, a^) defined by Eq. (4.34), and the body's 
position is evaluated on the past light cone v = constant so as to yield a meaningful 
comparison. Here also we find that the horizon tide (as measured by the displacement 
field) leads the source of the tide (as measured by the orbital position on the light 
cone) . 
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Figure 3. Shown in solid red is the polarization (t+ associated with the shear 
tensor of the horizon's intrinsic geometry, rescaled by a factor of m, as a function 
of i;/(2M). Shown in dashed blue is the polarization C+ associated with the Weyl 
tensor, rescaled by a factor of m/(2M), as a function of v/{2M). Both quantities 
are calculated at position o = ^ and /3 = on the horizon. The plots show clearly 
that the horizon tide (as measured by the shear tensor) leads the tidal field (as 
measured by the Weyl tensor) by a time interval of order = 4M. The Weyl 
tensor is noisy for early and late times because it is inaccurately computed by 
estimating the second derivative of 7_|_ with respect to v with finite-difference 
techniques. 



1.5 - 




-1 -0.5 0.5 1 1.5 2 

Figure 4. Snapshot of the tidal bulge at v/{2M) = 4.4874 as described by 
the dimensional tidal displacement p evaluated as a function of (5 on the black 
hole's equatorial plane a = ^. The figure shows, in the same x-y plane as in 
Fig. 1, the surface r = 2M of the unperturbed horizon (in thin black) as well 
as the surface r = 2M[1 + {M/m)p] (in thick red), which grossly exaggerates 
the horizon deformation by a factor of Mjm, to make it visible. The figure also 
shows (red disk) the position of the orbiting body at this value of advanced time 
v\ we have r/{2M) ~ 2.0273 and <j> ~ 0.6005, leading to the Cartesian positions 
x/(2M) ~ 1.6726 and y/{2M) ~ 1.1455. The tidal bulge and orbiting body are 
intersected by the same light cone v = constant, and here also we see the tidal 
bulge leading the source of the tidal field. 
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Appendix A. Late-time behaviour of horizon quantities 



Some of the horizon quantities (such as -y''™", ^odd^ a;''™", and w° ) can be expressed 
purely in terms of the current value of the master functions, while others (such as 
^trace^ ^even^ g^j^^ j^odd^ InvolvB integrals of the master functions. We wish to verify 
that all horizon quantities properly vanish at f = oo, assuming that ^'ovon(^', 2M) and 
\l/odd(^'5 2M) decay at least as fast as an inverse power law in v; this is the late-time 
behaviour expected of radiative tails that linger on after the external processes that 
produce the perturbation have shut down. 
The general structure of the integrals is 



x{v) = a;(t;o)e-'=°(^-'"') - f e'^-^^''-'''^ F{v') dv' 

Jvn 



(A.l) 



and for our purposes here we assume that the source fuuc;tion F{v) varies over a 
timescale r that is very long compared with Kq ^ = 4M. In the case of an inverse- 
power falloff, for example, we assume that vq is sufficiently large that F{v') oc (i"')"^ 
within the integral, with p > 0. Then F o: (w')"^"""^ and the timescale r can be 
identified with F/F oc v'; this is indeed much larger than AM for the specified domain 
of integration. In these circumstances we can evaluate the integral and express it as 
an asymptotic series in powers of {kot)~^ <C 1. If we let 



I[F] 



-Ko{v — v') 



F{v') dv', 



then the identity 



I[F] 



■I[F] 



(A.2) 



(A.3) 



follows immediately by integration by parts. Repeated applications yield 



1 r 



+ 



1 



I[F] = F{vo) - Ko'F{vo) + Ko^F{vo) + ■ 



-Koiv — Vo) 



F{v)-Ko'F{v) + Ko^F{v) + --- 



(A.4) 



in which each term within the square brackets is smaller than the preceding one by a 
factor of order (kqt)^^. With this we arrive at 

x{v) = — \kox{vo) + F{vo) - Ko^F{vo) + Ko^F{vo) + ■■■] e-''°(''-''°) 

Kq L J 

- — \f{v) - Ko^F{v) + KfF{v) + •••!. (A.5) 

/to L J 

At large v the first set of terms decay exponentially, and x{v) is dominated by the 
second set of terms. A good approximation is then x{v) ~ —Kq^F{v), and x is seen to 
decay at the same rate as F{v). This shows that our integrals are indeed well-behaved 
in the limit v ^ oo, and that the horizon quantities decay at the same rate as the 
master functions. 
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